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Vortex dynamics in inhomogeneous Bose-Einstein conden- 
sates are studied numerically in two and three dimensions. We 
simulate the precession of a single vortex around the center of 
a trapped condensate, and use the Magnus force to estimate 
the precession frequency. Vortex ring dynamics in a spherical 
trap are also simulated, and we discover that a ring undergoes 
oscillatory motion around a circle of maximum energy. The 
position of this locus is calculated as a function of the num- 
ber of condensed atoms. In the presence of dissipation, the 
amplitude of the oscillation will increase, eventually resulting 
in self-annihilation of the ring. 

PACS numbers: 03.75.Fi, 05.30. Jp, 67.40.Vs 



I. INTRODUCTION 

Among the most important phenomena associated 
with Bose-Einstein condensation (BEC) is the quantiza- 
tion of vorticity, which is intimately connected with the 
existence of persistent currents and superfluidity in quan- 
tum fluids. Study of quantized vortices has been confined 
mainly to liquid Hell where detailed comparison to 
mean-field theory is complicated by strong interactions 
between atoms. However, such considerations are much 
less important for the recently achieved BEC in atomic 
vapors H-^. In this case the condensate can be accu- 
rately described by the Gross-Pitaevskii (GP) equation, 
an example of a nonlinear Schrodinger equation, whose 
properties are well-known. This equation admits vortex 
solutions, where a non-zero circulation is accompanied by 
a zero in the condensate density. The density variation 
defines the vortex core, with a size of ~ 1 /zm (c/. ^ 1 A 
in Hell). Thus, vortices may be directly observable by ab- 
sorption imaging ||^,^ or other detection schemes [0-0 • 
Apart from their intrinsic interest, vortices play an im- 
portant role in the breakdown of superflow in Bose fluids 
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Despite the enormous progress of experiments on 
atomic condensates a vortex state has yet to be 

observed. One possible procedure involves rotating the 
condensate |lT^,15t at a rate exceeding a critical angular 
velocity |lq-^[lq , creating a singly-quantized vortex line 
along the axis of rotation [|l^,0. Larger angular veloc- 
ities might be expected to produce vortices with higher 
quanta of circulation. However, previous studies [ pT|j22| 
suggest that such vortices are unstable, leading to arrays 
of singly-charged vortices similar to those found 

in HeH fl. 



A separate but complementary idea is to use a tightly- 
focussed far-off resonant laser beam which creates 
vortex pairs when dragged through the condensate. Al- 
ternatively, by stirring the condensate angular momen- 
tum is transferred, and a single vortex may be produced 
pGf . A vortex ring may be formed by translating one 
condensate through another |^^, or by 3D soliton decay 

Experimentally, the Bose-condensed gas is usually con- 
fined in a magnetic trap, often modelled by a harmonic 
potential. This profoundly alters the condensate proper- 
ties in many ways [ p^ , pO[ , the most apparent of which is 
its spatial inhomogeneity. This results in a variation of 
vortex energy as a function of position, with a maximum 
at the center of a non-rotating trap. As a consequence, 
a single vortex precesses around the condensate center 
[^,01 • In addition, the vortex is thermodynamically un- 
stable, and dissipation at finite temperatures leads to its 
expulsion from the cloud |^,^. The instability is also 
apparent in the excitation spectrum, with the existence 
of a mode possessing riegative energy with respect to the 
ground state |^ , ^ , P5| . 

In this paper, we study the motion of vortices in 
trapped Bose condensates, by numerical solution of the 
GP equation. As discussed in Section II, this equation 
is valid in the limit of low temperatures, and describes 
the conservative motion of a vortex. Section HI presents 
measurements of the precession frequency of a single vor- 
tex in two and three dimensions, and compares the re- 
sults to analytical expressions. We use a Magnus force 
argument to estimate the precession frequency in 2D. In 
Section IV we study vortex rings, and find that they per- 
form a cyclical motion. A stationary state can be found 
which corresponds to a ring with maximum energy. This 
point of unstable equilibrium is equivalent to that of a 
single vortex in the center of a non-rotating condensate. 
Finally, we summarize in Section V, and briefly discuss 
finite temperature effects. 



II. THEORY 



A. The Gross-Pitaevskii equation 



In experiments on atomic vapours evaporative 
cooling can be extended to very low effective tempera- 
tures, such that the non-condensate fraction is very small. 
The densities are sufficiently low that interactions can be 
represented by a pseudopotential of the form UoS{r --r'), 
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where Uq = ATrh^a/m, and a is the s-wave scattering 
length. This leads to the Gross-Pitaevskii (GP) equation 
|p6|,p7t for the condensate wavefunction, \['(r,t), given by: 



ih 



'dt 



t2 

- — + y + 7V[/o|*|^ ) ^, 
im 



where N is the number of atoms, each of mass m. Conse- 
quently, A'^jvl/p is the condensate number density, which 
satisfies the normalization condition: 
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The harmonic trapping potential is denoted by ^(t') = 

It is convenient to scale in dimensionless units. We 
use harmonic oscillator units (h.o.u.) where the units 
of length, time and energy are (?i/2mwa;)^/^, uj~^ and hjjx 
respectively. We also consider a frame rotating about 
the z-axis with angular velocity 57, where the angular 
momentum operator is given by 
taining the normalization condition 



i{ydx -xdy). Re- 
), Eq. (P becomes: 
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where V = \{x^ + rjy'^ + ez^) and the anisotropy 
parameters are defined as 7y = uj^/uf^ and e — 
Lol/uj^. The interaction parameter is given by C = 
{NUo/huJx){2muJx/h)°'/^, where a is the number of di- 
mensions. In 2D, where atoms are confined in the x — y 
plane, TV represents the number of atoms per unit length 
along z. 

can be integrated by various numerical meth- 
As in our previous work [^5| , ^ , we utilize a 
Fast Fourier Transform (FFT) technique. Stationary so- 
lutions of (|^), can be represented by ^'(r,t) = \l/(r)e~*''* 
where fi is the chemical potential. They are found by 
propagating a trial wavefunction (e.g. a Gaussian so- 
lution for C = 0) in imaginary time, t — > —it. In 
contrast to real time, the resulting evolution operator 
is non-unitary, so the wavefunction must be renormal- 
ized after each time-step. The ratio of norms pro- 
vides a convenient estimate for the chemical potential, 
fi = (2Ai)-iln[(|*(i)|2)/(|^'(t Ai)|2)]. In imaginary 
time, excitations are exponentially damped, and both ^' 
and fi rapidly converge to a stationary solution, provid- 
ing accurate initial conditions for time-dependent simu- 
lations. 

The free energy per particle, E, is given by: 
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where for a stationary solution it can be shown that: 

E = t^-^j |*|*dV (5) 

Imaginary time propagation minimizes the chemical po- 
tential, /i. 



B. The Vortex State 



Using the Madelung transformation j^], the conden- 
sate wavefunction can be represented in terms of its den- 



(1) sity p(r,t) = |*(r,t)p and phase 5(r,t) by * = ^e 
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The superfluid velocity is given by Vg = fi.(^'*V^ — 
\I'V^*)/imp = {h/m)VS. So, the circulation around 
an arbitrary closed loop is: 
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where the property of a single-valued wavefunction re- 
stricts n to an integer value. From (^, V x. Vg — 
2'Kni5{r — roi), and the superfluid is irrotational ev- 
erywhere except for vortices at r = roi- We define a 
circulation vector k on the vortex axis, which has mag- 
nitude |k:| = K and direction V x Ds- 

A stationary state exists for a single vortex line at ro = 
0. If, for simplicity, we assume a non-rotating isotropic 
trap in the x — y plane, t hen this state can be represented 
as ^'(r, 0, z) — \/ p{r, z)e*"'^. Substitution into (Q) gives: 
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where f ~ {2mfj,)^^'^r/h, z 
asymptotic forms: 



= (2mpf/'^z/h, and / = 
Eq. (|^) yields the following 
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and p ~Q for r > i?^ and \z\ > Rz, where i?^ ~ 2p/niuj'^ 
and Rl = 2p/mu)'^. The parameter ^ — ft/(2m/i)^/^ = 
(SnpQa)^^^^ is the coherence length in the center of the 
condensate Jql , and determines the size of the vortex 
core. Eq. (§)is presented in [|7jl9|3^ and is analogous 



to the well-known Thomas-Fermi (TF) approximation for 
the ground state, valid for large NUq. In this limit, po 
corresponds to the density in the center of the ground 
state condensate. 

For the special case n = 1, Eqs. (|) and (|9|) can be 
interpolated by a simple function, whic : 
straight, off-axis vortex line at Tq is: 
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where p = when the right-hand side is negative. We 
compare Eq. ( [lo| ) to the numerical solution of the GP 
equation in Fig. 1. The latter is found by imaginary time 
propagation, imposing a phase variation of 2t: around the 
vortex line at each time-step. Eq. (|lO|) provides a good 
estimate especially in the high-A^ limit, in a similar way 
to the TF approximation to the ground state. 
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FIG. 1. Cross-section through a singly-quantized vortex 
line, showing condensate density as a function of position. 
The solid line plots the exact 3D wavefunction as calculated 
from imaginary time propagation of Eq. (^, for C — 4000, 
n = and e = rj — 1. The dashed line represents the ana- 
lytic approximation (0), with po ~ 2.443 x 10"^, ^ ~ 0.3199, 
ro = 0, and R± = Rz — 6.252. These parameters are cal- 
culated from the chemical potential ptf = (15C/647r)^''^, 
given by the normalization condition of the ground-state 
Thomas- Fermi wavefunction. 



III. SINGLE VORTEX MOTION 
A. Two Dimensions 

In this section we study the dynamics of a singly- 
quantized (n = 1) vortex line in a non-rotating conden- 
sate. Imaginary time propagation of the GP equation (|^) 
is used to provide an initial condition for the real-time 
simulation. If we consider a vortex at position ro relative 
to the axis of a non-rotating trap, then the free energy 
of the system, Eq. (jj), is found to attain a maximum 
when ro — (see Fig. ||). So, a vortex initially at the 
origin will remain stationary, but is unstable to infinites- 
imal displacements. The GP equation implies that the 
system is Hamiltonian, and therefore an off-axis vortex 
will follow a path of constant energy corresponding to 
precession around the trap center. The presence of dissi- 
pation will lead to drift towards lower energies, causing 
the vortex to spiral out of the condensate |3^,Q . 

If the trap is rotated with angular velocity, f2, then 
the energy of a central vortex, E^ot , decreases such that 



E^ot = En — n^ ||22| . Appearance of a vortex becomes en- 
ergetically favorable when E-^^t < Eq, so that the critical 
angular velocity is simply: 
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En — Eq 



In the TF limit 
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For \n\ > fie, an on-axis vortex attains global stability; 
however, there remains an energy barrier for vortices en- 
tering from the edge (see Fig. ||) . 

Inspection of Fig. Halso reveals that that above an 
angular velocity flmi the vortex attains a local minimum, 
so a vortex is metastable when < \n\ < ftc- In the TF 
limit, f2„i = 3f2c/5 |31|, while in the non-interacting limit 
flm |^l|- In both the weak and strong coupling 

limits, flm — — Wa, where uja is the frequency of the so- 
called 'anomalous' mode obtained from solution of the 



Bogoliubov equations |31 3^,^,^. It is thought that to a 
corresponds to the frequency of vortex precession, thus 
linking vortex dynamics and instability. 
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FIG. 2. Energy, E, as a function of n = 1 vortex displace- 
ment, in a 2D condensate rotating with angular frequency, 
Q, {C = 1000, ?7 = 1). The top solid curve corresponds to 

= 0, where increases in steps of 0.05 as one moves to- 
wards the lowest curve. The dashed line marks the energy of 
the condensate without a vortex. 

First, we consider the simplest case of a vortex in 
two dimensions, which corresponds experimentally to a 
condensate confined in an axisymmetric cylindrical trap 
(where e ^ 0, 77 = 1). For = 0, simulations show 
that an off-center vortex accelerates from its initial con- 
dition, soon attaining a near-constant angular velocity cu 
around the trap center, such that the instantaneous ve- 
locity is t>L = ijJK. X r. The angular velocity, iv, is plotted 
as a function of interaction strength and initial position 
in Fig. ||. For small C, uj is averaged over a few cycles 
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(e.g. 3 revolutions for C = 200). However, for higher 
C, numerical instabilities can restrict the simulations to 
less than a half-cycle: the error bars reflect the resulting 
uncertainty. 




100 



1000 



10000 



FIG. 3. Vortex precession frequency, uj, in a 2D condensate 
at ro = 0.5, as a function of interaction parameter C = SnNa. 
Filled circles show the results of numerical simulations. The 
triangles indicate the Magnus force estimates, obtained from 
the gradient of the numerical values of Eq. The analytical 
Magnus force estimate, Eq. (|l5|), and Eq. ( |l^ ) are plotted with 
dot-dashed and dashed lines, respectively. The TF vortex 
metastability frequency, | Ai5, is plotted as a solid line. Inset: 
(i) as a function of vortex position ro, for C = 1000. 

The precession frequency decreases as a function of 
increasing interaction strength, C. An intuitive semi- 
analytical argument for this behavior can be formulated 
in terms of the Magnus effect, familiar from classical 
hydrodynamics, and in superfluids and superconductors 
When the background fluid flows past the cir- 
culating fluid connected with the vortex, a pressure im- 
balance is created perpendicular to the direction of the 
background flow. The resulting Magnus force must bal- 
ance the force due to the variation of energy E with posi- 
tion, i.e. dE/dro = mpK xvl, where vl is the velocity of 
the vortex line relative to the ambient condensate. Note 
that the Magnus force can also be produced by flow of the 
condensate around a stationary vortex, as in a rotating 
trap. So, one expects that lo ~ n„i- 

We find E by evaluating the functional, Eq. with 
a wavefunction grown in imaginary time with 17 = 0. 
Using the Madelung transformation, the first term splits 
into a 'quantum pressure' and a 'kinetic energy' term: 
iV^fp = (^y/pf + p(VS')2. A numerical differentiation 
of E with respect to tq gives an estimate for u using the 
Magnus force argument (plotted as triangles in Fig. ||). 
However, these estimates are sensitive to small numerical 
errors in the energy. 

To obtain an analytical estimate, we observe that /i is 
approximately constant at small rg and high C. This 



may be shown by using the decomposition ^{x,y) = 
$(x, ?/)G(x, y), where an off-set vortex core ($) is im- 
printed on the TF ground-state {&), as implied by Eq. 
(|l0|). Starting from the GP equation (|^), and noting that 
in the TF limit is slowly varying, then it is simple to 
show that: 
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where /itf is the chemical potential of O. As the Lapla- 
cian is spatially invariant, it follows that fi is independent 
of the offset at small rg. This can also be justified nu- 
merically, though the approximation only becomes valid 
at high C. Using Eq. (0) it follows that: 
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Substituting Eq. ([l^) for p then gives an estimate for the 
precession frequency (to logarithmic accuracy) : 
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This result may be compared to the expression ob- 
tained by Svidzinsky and Fetter pl[ |, using a time- 
dependent variational analysis: 
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again valid for small tq . These expressions are plotted to- 
gether with the numerical results in Fig. ^. In addition, 
we plot U,jn = where IS.E is the energy difference 

between the ground and n = 1 vortex states. Recall that 
one expects that oj ~ in the TF limit. The numerical 
results lie between the results of (|l^ ) and the metastabil- 
ity curve. All of the curves reproduce the correct func- 
tional dependence at high C. Note that these expressions 
are only valid for small ro; the vortex precesses faster as 
it nears the edge of the condensate, as shown in Fig. ^ 
(inset). 

Compressibility effects become important when the 
vortex is accelerating or when the velocity is an apprecia- 
ble fraction of the speed of the sound, Cg = {pNUo/my^'^. 
In an infinite compressible fluid, phonons may be emitted 
by a moving vortex, leading to a drift of the vortex to 
lower energies However in a finite condensate, where 
excitations remain confined in the region of the vortex, 
no net drift is expected. At the beginning of the sim- 
ulations, we observe an increase in radius of precession, 
together with excitation of an elliptical center of mass 
mode at the trap frequency. In addition, surface waves 
are created when the vortex is near the condensate edge. 
However, as expected we do not observe a sustained vor- 
tex drift (for T up to 110, corresponding to ~ 6 full 
cycles for C = 200). A drift to lower energies would be 
expected where a thermal cloud damps the motion (i.e. 
at finite temperatures). Nevertheless, the vortex lifetime 
is expected to be long, especially for large numbers of 
atoms 1391. 
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B. Three Dimensions 

Vortex dynamics become more complex in 3D, as the 
vortex line can deform along its length. In classical and 
quantum fluids, this can result in the propagation of he- 
lical waves along the line — so-called Kelvin modes . 
In simulations of 3D vortex motion, we have observed 
line deformation and oscillations. However, the inhomo- 
geneity of the condensate complicates matters, and the 
motion is difficult to resolve into simple Kelvin waves 
characteristic of the bulk condensate. The amplitude of 
the oscillations are typically small, and as a consequence 
helical waves are likely to be difficult to detect experi- 
mentally. Moreover, it is worth noting that the energy of 
the vortex increases as it lengthens. Hence, in the pres- 
ence of the dissipation the line will tend to straighten, 
effectively damping the Kelvin modes. 

Fig. ^ compares numerically measured values of the 
precessional frequency with the TF result of Svidzinsky 
and Fetter Eq. (|l6l) . It can be seen that the fre- 
quency dependence is well described; however, the nu- 
merical results are significantly higher 20%), but con- 
verge slowly to the analytical expression towards high C. 
This disparity may be due to effects resulting from the 
curvature of the line, which will modify predictions that 
assume a rigid line motion. Numerical values of the TF 
metastability frequency, |Ai?, are also found to be lower 
than the observed precession frequencies. 
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FIG. 4. Vortex line frequency, u, in a 3D condensate 
(77 = 1, e = 9) at ro — 0.5, plotted as a function of inter- 
action parameter, C = 8TvNa{2muj/h)^^^ . Squares display 
numerical results, while the dashed line plots Eq. (^6|). The 
solid line shows the TF metastability frequency, |Ai5. 

IV. VORTEX RING MOTION 

The motion of a vortex ring in a trapped BEC may 
be understood in terms of a sum of two contributions 
to the velocity of each element in the ring. First, the 



precession due to the inhomogeneity of the condensate, 
as discussed for a single vortex in Sec. HI, and second, the 
velocity induced by the remainder of the ring, Wi„, which 
is directed along its axis (defined as the z-axis). For a 
spherical condensate, the total velocity on each element 
is given by: 

V = VinZ + t^k X r, (17) 

where k, defines the direction of the circulation at the 
element, and lo is the precession frequency. In a homoge- 
neous Bose fluid i^in = (ft/2TOi?^)[ln(8i?,./C) - 0.615] @, 
where Rr is the ring radius. Consider a ring at z = 0, 
r = Rr- If the radius is small, the induced velocity dom- 
inates and the ring moves in the +2 direction, while if 
Rr is large the precession dominates and it travels back- 
wards. In addition, the precessional term leads to ring 
expansion for z > and contraction for z < 0. Thus, the 
two terms produce an oscillatory motion of the ring. 
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FIG. 5. Vortex ring energy in a cylindrically symmetric 
condensate (77 = e = 1, C = 2000) as a function of radius, 
r, and z-position. The energy contours are equally spaced 
between 5.6179 and 6.1159. The bold line shows the motion 
of one element of the vortex ring, where the circles represent 
the position at equally-spaced times (every T = 1). The ring 
begins at (1.5, 0), marked by an arrow, and cycles around the 
energy maximum in a clockwise direction. 

One can also understand the ring motion as a trajec- 
tory around an energy maximum, in analogy with the 
single line vortex. To demonstrate this, in Fig. [s] we plot 
the energy of an on-axis ring as a function of its radius, 
r, and z-position. Without dissipation one would expect 
the ring motion to follow an energy contour; however as 
is apparent in Fig. ^ this is not exactly true. Accelera- 
tion of the ring at the beginning of its motion results in a 
back-action on the condensate, exciting a center-of-mass 
mode, and the subsequent ring dynamics are complicated 
by the underlying motion of the condensate. In addition, 
for small C we observe a decay of the ring to lower en- 
ergies over the first cycle of its motion. As for a single 
vortex, this effect is associated with the compressibility. 
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which results in acoustic emission from the moving ring 

The energy maximum at r = Rcq, z = 0, corresponds 
to the point where the two velocity contributions in (|l^ ) 
are equal and opposite, leading to a ring in unstable equi- 
librium. To obtain an analytical estimate for this posi- 
tion, we approximate the ring energy by taking the en- 
ergy of a single 2D vortex, and integrating around a 
circle of radius Rr, such that Er = 2TrRrEy. The domi- 
nant contribution to Ey is given by the kinetic energy, so 
Ey ~ (™/2) / pWgd'^r. Taking Vg = k x r/27rr^, where 
we translate the cylindrical coordinate system so that the 
origin lies on the vortex axis, and using (QOh gives: 



Er 



2n-^h^P()Rr 
m 



Rl 



(18) 



neglecting terms of order and higher. For a ring of 
radius Ry at Zo, Tq = Rl + z^. This expression describes 
the qualitative features of Fig. ^; however, it tends to 
over-estimate the energy near to the peak 10% at 
C = 2000) and is a poor approximation as Ry 0. Nev- 
ertheless, we can obtain an estimate for the equilibrium 
position from Eq. (n8|), which yields zo = and: 
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where /3 = (muj/h)^. In the TF hmit (_Rj_ — > oo), R^q ~ 
R±/V3 ~ 0.577i?^, which is close to the results from 
numerical solution, where Rcq ^ 0.54i?_L at high C (Fig. 
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FIG. 6. The equilibrium ring radius, i?eq, as a func- 
tion of nonlinear parameter, C = 8nNa(2mu>/hy^'^ , in 
a spherically symmetric condensate. The lower plot 
shows the ratio of R^q to the Thomas-Fermi radius, 
Rtf = R± = {2^TF /mui^y''^ , where the chemical potential is 
given by the normalization condition of the TF wavefunction, 
such that fiTF = {15C/6Anf^^. 

An experimental technique for ring production was 
proposed in Ref. The method utilizes a two- 

component BEC, such that when the smaller component, 
|2), is translated with respect to the other, vortex rings 
are created in the larger condensate |1). 

To model the two condensates, a pair of coupled GP 
equations are solved for the wavefunctions \l/i(r,i) and 

idt^, = [-V^ + V, + C + iM/^f )] (20) 

where i,j = 1,2 (i ^ j) |^^. The initial states of 
the simulations are created by imaginary time propaga- 
tion as previously, where the normalizations are set so 
that a variable fraction of atoms are in each condensate, 
X = (|4'iP)/(|4'2p), and (|*i|2)-t-(|*2|2) = 1. Our simu- 
lations then follow the creation and subsequent dynamics 
of the vortex ring, an example of which is shown in Fig. 
0. The trajectories roughly follow a contour of constant 
energy (see Fig. I). The ring is created from zero radius 
at t w 1.6 and y « 1.8. It then expands and travels 
forward, before turning and progressing backwards along 
the edge of the condensate. Finally, it turns again and 
collapses to a point, where the ring is annihilated. The 
annihilation produces a sound wave, which decreases in 
amplitude as it propagates along the z-direction. The 
sound wave then disappears at the edge of the cloud. 

In the presence of dissipation, annihilation will even- 
tually occur for any ring as a culmination of a decay to 
lower energies. This is equivalent to the instability mech- 
anism for a single vortex line (Sec. III). 

Vortex ring detection is likely to present considerable 
challenges to experimentalists. The simplest and most 
widely used method of condensate imaging is by mea- 
surement of probe laser absorption, after release of the 
condensate from the trap. The subsequent ballistic ex- 
pansion results in an effective magnification of a vortex 
line, which should then be visible as a density 'hole' [1). 
However, for rings the density minima are obscured by 
the rest of the condensate along any line-of-sight. One 
solution is to view slices of the condensate after ballistic 
expansion, using light sheets. An alternative method is to 
study the center-of-mass motion of coupled condensates, 
yielding details of the mutual drag that reveal vortex ring 
formation Collective excitations are also utilized in 
another scheme, proposed in [^,0, where a vortex line 
splits the degeneracy of the quadrupole mode of the con- 
densate. However, further work is needed to extend this 
analysis to vortex rings. 
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FIG. 7. Vortex ring motion in a condensate (77 = e = 1) 
after creation from an object. The trajectory of the ring is 
determined by solving Eq. @, with v = 1.75, C = 1100 
and X ~ 10/11. The upper plot shows the ring radius (solid 
line) and z-coordinate (dashed) as a function of time, while a 
parametric plot (bottom) displays the ring radius, r, against 
position z. 



V. SUMMARY 

We have studied the motion of vortex Hnes and rings in 
Bose-Einstein condensates in harmonic traps, by numeri- 
cal solution of the Gross-Pitaevskii equation. We consid- 
ered a single vortex in two and three dimensions. At the 
center of a non-rotating condensate the vortex state pos- 
sesses maximum energy, corresponding to unstable equi- 
librium, while an off-center vortex undergoes precession 
around this maximum. The precession frequency was 
measured and compared to theoretical models. 

Vortex rings were also found to undergo an oscillatory 
motion. For a particular radius, the ring energy is a max- 
imum, corresponding to state of unstable equilibrium. A 
lower energy ring precesses around the locus of maximum 
energy. Also, a ring which is created at a point will even- 
tually collapse to a point, resulting in self-annihilation. 

The study of dissipative vortex dynamics at finite tem- 
peratures is particularly interesting. In this case, inter- 
action of the vortex with the thermal cloud leads to a 
transfer of energy, and consequently a decay of the vor- 
tex state. One can consider this to be due to anisotropic 
scattering of thermal excitations , where the momen- 
tum transfer is apparent as a frictional force on the vor- 
tex. The transverse component of this force results in 



expulsion of the vortex from the condensate. A full mi- 
croscopic model of vortex dynamics is needed to describe 
this behaviour and will form the focus of future work. 
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